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ABSTRACT 

An analysis is made of the effect of sea state on Doppler-radar 
measurement of aircraft ground speed. Previously published data taken 
at both microwave and visible-light wavelengths are used to estimate 
the ground-speed error in a representative situation. The magnitude of 
the error is determined by the sea state and by the flight direction 
relative to the ground-wind direction. The results indicate that a 
probable error of less than one percent is generally achievable after 
correction for the effect of sea state. 
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SUMMARY 

An analysis has been made of the effect of sea state on Doppler-radar 
measurement of aircraft ground speed. Previously published data taken at 
both microwave and visible-light wavelengths were used to estimate the 
ground-speed error in a representative situation. The magnitude of the 
error is determined by the sea state and by the flight direction relative 
to the ground-wind direction. 

The results indicate that a probable error of less than one percent 
is generally achievable if one estimates the most probable value of the 
error for each heading, and for the given sea state, and then applies 
this estimate as a correction. 

INTRODUCTION 

In some meteorological experiments covering a large area over the 
ocean, it may be necessary to measure the wind speed and direction at 
several altitudes. This can be done, using instrumented aircraft, by 
combining Doppler radar measurement of ground speed with the aircraft's 
airspeed. The computed wind speed is the difference of these two much 
larger numbers, both of which are subject to error. The effect of sea 
state on the error in ground speed measurement by Doppler radar will be 
treated here. The analysis relies on previously-published experimental 
data taken at both microwave and visible-light wavelengths. 

This analysis was initially performed to assist in the planning for 
project BOMEX (Barbados Oceanographic and Meteorological Experiment), 
which was conducted in May-July 1 969 by ESSA with the participation and 
assistance of numerous governmental and university organizations. 
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Doppler radar measurement of aircraft ground speed over the sea is 

•5f 

subject to two errors not found over land. Both are due to the character 
of the sea surface and vary quantitatively with the sea state. The first 
error is due to the fact that the radar beam is reflected by the sea 
surface in a different manner than that in which it is reflected by land. 
The error, which may be termed a "calibration error," may be on the order 
of one percent of the indicated ground speed. The second error is due 
to the fact that the radar-reflecting sea surface may, under some circum- 
stances, appear to have a velocity, relative to "true ground," whose 
magnitude depends on the surface wind speed. This error, which may be 
termed a "surface-speed error," is some fraction of the surface-wind speed. 
Both errors depend on the angle between the aircraft heading and the wind 
direction at the sea surface. 

The dependence of the calibration error on the nature of the sea 
state appears fairly well documented by several observers (references 1-4) 
who have reported actual Doppler-radar data. Unfortunately, these data 
are given without specifying the angle between the radar beam and the 
surface-wind direction. The dependence of the calibration error on the 
angle between aircraft heading and surface wind ("crosswind" vs "upwind- 
downwind" ) is less fully documented. We have had to use the data reported 
by references 5 and 6, taken at visible-light wavelengths. There is no 
evidence that these data are inapplicable to radar wavelengths; however, 
the data cover a narrower range of sea states than was covered by 
references 1-4. 

The dependence of the surface-speed error on the surface-wind speed 
is presented by references 1, 2, and 7* However, it is not certain that 
these represent independent determinations; they may partially stem from 
the same source. Furthermore, the magnitude of the surface-speed error 
is stated in a less precise form than is the calibration error. We have 
had to apply some intuition in order to make the best possible estimate 
of the surface-speed error. 

* The simple error due to ocean current (l to 2 knots) is neglected in 
this analysis. 
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In the following analysis, we first discuss the measurement of the 
velocity of the aircraft relative to the sea surface, and derive the 
corresponding calibration error. Then, we discuss the determination of 
the velocity of the sea surface relative to "true ground" and derive 
the corresponding surface-speed error. 

A numerical example illustrates how these errors are combined. 


System Configuration 


For simplicity, the system to be considered will be the Janus 
system (reference 9 ) shown in Fig. 1. The simplification is tantamount 
to assuming zero drift angle, so that aircraft heading coincides with 
the track of the aircraft over the ground, and to assuming that the plane 
through the forward and rearward beams includes this ground track. The 
aircraft heading is in the x-direction. The horizontal plane through the 
mean sea surface is the x-y plane. Two radar beams from the antenna at 
(0, 0, h) reach the sea surface and are reflected back to the antenna. 

Each beam is very narrow, rectangular in cross section, with included 
angle 2a in the x-z plane and included angle 2(3 in the orthogonal direction. 
The beams intersect the x-y plane in the two rectangles shown. The forward 
and rearward beams both make a mean angle Q q with the z-axis in the x-z 
plane. If (3 lies along an isodop, the beam can be treated as two- 
dimensional, so that only what happens in the x-z plane need be considered. 

Infinitely narrow beam. For the idealized condition 0=0, the basic 
Doppler formula holds 


Af 

f 


4v 

— sxn 0 
c o 


( 1 ) 


where Af/f is the fractional frequency shift, c is the speed of light, 
and v is the component of aircraft velocity relative to the sea surface 
in the direction of the aircraft heading. No frequency shift is produced 
by relative-velocity components in the y- or z -direction. 

Beam of nonzero width. For the realistic case a > 0, the_ sin 0 q in 
Eq. (l) must be replaced by some average over the interval 

(e - a) < e < (e Q + a). 


( 2 ) 
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Assume that the antenna beam has a uniform radiated-power distribution 
over this interval, but that the power reflected at angle 0, per unit 
element of area of the sea surface, varies with the angle of incidence 0. 
Then, if cr°(g) is the back-scattering cross section of the sea surface 
per unit area of sea surface, Eq. (l) is replaced by 


Af 

f 


4v 

c 


/ 


cr°(0 ) sin 0 


o°(0) • d0 


- d0 


(3) 


This formula assumes that cr°(0) is the same for both forward and 
rearward beams. For the situations treated in this analysis, this 
assumption appears to be valid (references 5 and 6). 

Distinction Between Land and Sea 

Over land covered with reasonably dry vegetation, the value of o°(e) 
is independent of 0 (references 2 and 4), so that the ratio of the 
integrals, that appears in Eq_. (3), is very close to sin 0 . However, 
over the sea, one encounters two types of radar reflection: 

Specular reflection from smooth waves. Here, the statistical 
distribution of wave slopes is approximately Gaussian, and cr°(0) is a 
rapidly- changing function of 0 (references 1-4). The horizontal wave- 
propagation velocity has negligible Doppler effect in the Janus system 
because the water displacement is effectively a vertical oscillation 
(reference l). 

Scattering by whitecaps and spray. Whitecaps appear at wind speeds 
above 7 knots (Beaufort No. 3)j spray appears at wind speeds above 17 knots 
(Beaufort No. 5) (reference 8). The resultant scattering produces a value 
of a°(0) that is effectively constant over the interval Eq_. (2), just as 
in flight over land. However, the horizontal velocity of whitecaps and 
spray (relative to the "ground") causes a proportional Doppler frequency 
shift; spray, being wind-borne, has a higher velocity than whitecaps. The 
component of this horizontal velocity in the direction of the aircraft 
heading represents the "surface-speed error." 
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Combination of "both types of reflection. Both types of reflection 
can be included in Eq. (3) if we write 

a°(e) = a°(e) + a°(0) (*0 

where subscript 1 represents reflection by smooth waves and subscript 2 
represents reflection by whitecaps and spray. Note that 0°(0) is a 
constant. 

Error Due to Reflection by Sea 

Comparison of Eq„ (3) for sea reflection with the corresponding 
formula for land reflection (where 0° is a constant) leads to the 
fractional error in the fractional Doppler frequency shift. 



References 1 -b present radar data on a°(e) for various values of 

0 and various sea states. For the interval, Eq. (2), in which Oi « 0 , 
o 0 

the value of cr° can be approximated by an equation of the form 

0 0 = A + B0 (6) 

where B = (d0°/d0) 

C7—C7 


Substitution of Eq„ (6) into Eq. (5) yields 
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References 1 and 7 indicate that this calibration error can be partially 
compensated by a land/sea switch on a radar set. Such compensation 
reduces the probable calibration error by a factor of 2 (or more) through 
the expedient of inserting a correction equal to the most probable value 
of the calibration error , corresponding to the most probable sea state. 

Error Due to Surface Wind 

The calibration error just discussed affects the determination of 

the velocity v g of the airplane relative to the sea surface. There 

remains to be determined the velocity v of the sea surface relative 

s 

to "ground.” This latter velocity can be estimated from two data: 

(1) observation of the sea state 

( 2 ) independently-measured or estimated surface wind velocity. 

If (l) shows the sea to have only smooth waves (few or no whitecaps, 

Beaufort Wo. less than 4), the velocity v of the sea relative to 

s 

"ground" may be assumed to be zero. 

If (l) shows the sea to have numerous whitecaps (Beaufort Wo. 4 or 

higher), the velocity v of the sea relative to "ground" may be taken 
1 I s 

to be between and — of the surface wind velocity, according to 
references 1, 2, and 7- Presumably the fraction will be larger at higher 
surface-wind speeds. 

The component of v in the direction of aircraft heading, when 

s 

added to v , yields the aircraft ground speed in the direction of aircraft 
3 

heading. 

Calculations of Errors 

Radar scattering cross section of the sea surface. Graphs of cr° 
against sea state are shown for X-band (3 cm) waves by references 1 and 2 
in terms of Beaufort number, and by references 3 and 4 in terms of wind 
speed. Measurements of reference 4 were on a river instead of open sea. 
Angle between wind direction at the water surface and the radar-beam or 
aircraft heading is not stated. 

To indicate the order of magnitude of errors that might be encountered, 
computations will be made for the case 

0 = 20 °; 0 = 1 ° 
o 
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A linear approximation of the form a° = A + BQ to the published 
graphs , at Q = 20 degrees , gives the values of 



that are listed in Table I. The fractional error S(Af/f)/(Af/f ) for 
d = 1 degree , computed by Eq. (8), is also listed in Table I and is 
plotted in Fig. 2. 

Visible-light reflection by the sea surface. Because of the sparse- 
ness of data on cr° at radar wavelengths, it is useful also to use 

optical measurements of wave-surface slope and to derive an equivalent 
1 df 0 

value of — — • -r— therefrom. It is reasonable to assume that radar 

o° be 

waves and light waves are reflected in similar manner, since the wave 
length of either is very short compared to the wavelengths of the sea waves. 

References 5 and 6 report on specular reflection of light waves from 
the sea surface. They show that the statistical distribution of sea-wave 
slopes is Gaussian to a first approximation, with a smaller variance in 
the crosswind direction than in the upwind-downwind direction (Fig. 3 )* 

If the probability dP that wave slope lies in the interval dt]/ is 
given by 


dP = p(t) ' dij/ 

where 


( 9 ) 


p(t) = 7= 

2 it 



( 10 ) 


_2 

and \|/ is the variance of the distribution, then the reflectance R(\]/ q ) 
of the surface at a mean angle t 0 is proportional to the probability 


p(\|/ Q ) • At 


( 11 ) 


where Ai|r is the very small angle subtended by the transmitter-receiver 

(antenna) at the sea surface. A linear approximation to the curve of 

R(t) versus tb in the vicinity of t Q .> and i n the narrow interval between 

t + and t - , yields 

w o 2 Y o 2 ’ 



RU 0 ) = c i ° p ^o } * 



^ ' W t=t 
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( 12 ) 

( 13 ) 

( 14 ) 


where c. is a constant of proportionality that disappears from the final 
result (l 4 ). 

The quantity y • (~j can replace the quantity 


1 


0 ( 0 O ) 


(i-1 


e=e 


in Eq. (8), if 0 q = \)r « 


Figure 2 includes the errors so computed, using values of \|r from 
Fig* 3. The range of wind speeds covered is the range over which 
references 5 and 6 consider the data reliable. The maximum wind speed is 
15 knots, where whitecaps become numerous; the minimum wind speed, imposed 
by experimental limitations, is that where f becomes less than 0.4 i|r Q . 

Comparison between crosswind and upwind -downwind directions. The 
only reported data, presented in Fig. 2, that unequivocally indicates a 
distinction between the crosswind and the upwind-downwind directions of 
flight, are the light-wave data of references 5 and 6 at 15 knots. The 
error in crosswind flight is about 1.5 times the error in upwind-downwind 
flight . 

Based on this reported data, we make the following intuitive 
assumptions concerning radar waves: 

(1) At wind speeds higher than 15 knots, where whitecaps and spray 
provide stronger reflections, the ratio between crosswind and 
upwind-downwind directions should be closer to unity. 

(2) The radar data of references 1 , 2 , 3, and 4 were probably taken 
in the upwind-downwind direction. 

Numerical example. Table II presents numerical calculations of the 


errors under the following assumptions: 
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(1) 

( 2 ) 
(3) 

W 


(5) 


Aircraft speed is 300 knots 

0 = 20°, a = 1° 
o 

Data of references 1 and 3 are used for the upwind-downwind 
calibration error. 

Crosswind calibration error is assumed 1.5 times the upwind- 
downwind calibration error at wind speeds up to 15 knots; the 
multiplier is assumed to be 1.25 at 2 6 knots. 


Magnitude of surface-speed error in direction of the surface 
wind is assumed to be g of surface-wind speed at 15 knots and 
to be -jr of surface-wind speed at 2 6 knots. 

Note that the calibration errors would vary proportionately if 
aircraft speed were other than 300 knots, but that the surface-speed 
errors would not change. 

SUMMARY OF RESULTS 


The error in ground-speed measurement is the sum of a calibration 
error and a surface-speed error. The absolute magnitude of the calibration 
error at 300 knots is listed in Table II; this absolute error would vary 
in proportion to the speed. The surface-speed error, as listed in Table 
II, would not vary with speed. 

Analysis of the data in Table II then shows that for the conditions 
0^ = 20°, a = 1°, airspeed 300 knots, a probable error of less than one 
percent of the speed would be achieved if one estimated the most probable 
value of the ground-speed error from Table II and then applied this 
estimate as a correction. 
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Calibration Error Over Sea Surface, 

Relative to Land. Surface Calibration, 
A o „o 


Eor 9 = 20 


a 


1 


Reference 

Beaufort 

lumber 

Surface-wind 
speed, knots 

1 • K) 

a °(g o ) \&e J e=e ° 

Degree ^ 


1 

1 

1-3 

-o. 4 o 

-o.oo 64 


2 

4-6 

-0.30 

-o.oo 48 


3 

7-10 

-0.26 

- 0.0042 


4 

11-16 

- 0.24 

-0.0039 

2 

1 

1-3 

-0.32 

-0.0051 


2-3 

7 

-0o 22 

-0.0035 

3 


16 

-0,22 

-0.0035 


— 

26 

-0.13 

-0.0021 

4 


0-5 

-0.30 

-o.oo 48 


-- 

5-10 

- 0.40 

-o.oo 64 


-- 

10-15 

-0.20 

-0.0032 



15-20 

- 0.23 

-0.0037 
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Estimated Ground Speed Errors 

6 = 20 °, a = 1 ° 

o 7 


Beaufort 

Sea 

Surface-wind 

Type of 

Error in 

stated flight direction, knots 

Wo. 

State 

speed, knots 

error 

Upwind 

Downwind 

Crosswind 

1 

smooth 

2 

Calibration^ a ' 

-1.9 

-1.9 

-2.9 


waves 


Surface-speed 

0 

0 

0 




Total 

-1*9 

-1.9 

-2.9 

4 

whitecaps 

15 

Calibration^ a; 

- 1 * 1 * 

-1-1* 

-l.T 




Surface-speed 

1.9 

-1.9 

0 




Total 

0.8 

-3-0 

-l.T 

6 

spray 

26 

(a' 

Calibration' ' 

-0.6 

-0.6 

-5.2; 

* 

-0.75 




Surface-speed 

5-2* 

0 * 




Total 

4.6 

-5.8 

-0.75 


Two significant figures are not warranted. 

Q 

At 300 knots aircraft speed. 

















Fractional error 

5(Af/f)/(Af/f) 


RADAR OPTICAL 

• Ref. 1 A Ref. 5, upwind-downwind 

■ 2 < 5, crosswind 

k 3 V 6, upwind-downwind 

+ 4 > 6, crosswind 

Wind speed, meters per second 
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Wind speed, knots 


Figure 2. Error of Doppler frequency shift. 6 Q = 20°, a *= 1° 
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Figure 3. Optically- measured wave slope dispersion 


